Quantum correlations which imply causation 
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In ordinary, non-relativistic, quantum physics, time enters only as a parameter and not as an observable (T|: 
a state of a physical system is specified at a given time and then evolved according to the prescribed dynamics. 
While the state can, and usually does, extend across all space, it is only defined at one instant of time, in conflict 
with special relativity where space and time are treated on an equal footing. Here we ask what would happen 
if we defined the notion of the quantum density matrix for multiple spatial and temporal measurements. We 
introduce the concept of a pseudo-density matrix which treats space and time indiscriminately. This matrix in 
general fails to be positive for timelike separated measurements, motivating us to define a measure of causality 
that discriminates between spacelike and timelike correlations. Important properties of this measure, such as 
monotonicity under local operations, are proved. Two qubit NMR experiments are presented that illustrate how 
a temporal pseudo-density matrix approaches a genuinely allowed density matrix as the amount of decoherence 
is increased between two consecutive measurements. 



Ever since the pioneering work of Bell [2], the study of 
quantum correlations has proved fertile ground for gaining in- 
sight into fundamental physics. Much of that progress has 
been focused on spatial correlations, in the form of entan- 
glement and quantum discord [3-5], but a number of authors 
have extended this approach into the time domain. In partic- 
ular Leggett and Garg showed that quantum systems exhibit 
a form of timelike correlation which cannot be accounted for 
by any macro-realistic theory |6|. Similarly it has been shown 
that assumptions of realism and locality in time lead to a form 
of temporal Bell inequality, which again can be violated by 
quantum systems (7). Here we take a different approach: as- 
suming quantum mechanics a priori, and examining the cor- 
relations which can arise. Quantum states which violate the 
Leggett-Garg inequality necessarily exhibit the causal corre- 
lations we identify, and hence recent experimental demonstra- 
tions of violations of such inequalities may constitute a limited 
observation of such causal correlations I18 W21 . 

In quantum mechanics, each system in a multi-system 
quantum state is assigned a separate Hilbert space and these 
spaces are connected through the tensor product structure. 
The tensor product indicates that these systems are to be 
treated separately, though the joint state is also well defined 
at each instant of time. Here we explore extending this notion 
to different instances in time and assigning a Hilbert space 
to each different instant in time in much the same way as it 
is done in space. The resulting spatio-temporal state is then 
investigated. 

First we introduce the standard density matrix in quantum 
physics for qubits, although our ideas apply to subsystems of 
any dimensionality. We then show how to extend the con- 
cept of the spatial density matrix to different instances in time. 
The difference between spatial and temporal correlations is 
investigated through the introduction of the causality mono- 
tone, which is meant to capture the degree of "temporalness" 



in any quantum correlations. Finally, we present experiments 
using an NMR implementation that illustrate the basic prop- 
erties of the pseudo-density matrices. Our result shows that 
the simple phrase "correlation does not imply causation", fre- 
quently heard among philosophers of science, is incorrect if 
taken at face value. Our proposal to treat spatial and temporal 
correlations within the same quantum formalism clearly still 
discriminates between the two, albeit imperfectly: when the 
pseudo-density matrix fails to be positive, this means that it 
necessarily contains a temporal element; the converse of this 
is not true, as the pseudo-density can be positive without im- 
plying spacelike separation. 

The density matrix can be viewed as a probability distribu- 
tion over pure states, with p = X^PilV'iXV'il' where pi is the 
probability of the pure state \ipi) occurring. Given a density 
matrix p, the expectation value of a particular Pauli operator 
P is (P) = Ti(Pp). As the n qubit Pauli operators along with 
the identity form a basis for the space of Hermitian operators, 
and any density matrix p is necessarily Hermitian, it follows 
that any p can be written as p = a a I + J^. ctiPi, where Pi 

is the ith Pauli operator on n qubits, and {ai}™ = o me rea ^ 
numbers. Further, since Pauli operators are traceless, and all 
density matrices have unit trace, we have arj = 1/2™ and the 
expectation value for Pj is then given by 

(Pi) = Tr (Pj f^+E a ^ ) = Ta r C 1 ) 

Thus we have an alternate formulation of the density matrix 
in terms of the expectation value of Pauli operators, 

pA+Y^^p, ( 2) 

t=i 
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As we are interested in discussing correlations, we can express 
each n qubit Pauli operator as the product of single qubit op- 
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erators, yielding 
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where a = I, ai = X, er 2 = Y and cr 3 = Z. 

This above equation can be taken as definig a generalization 
of the density matrix. We consider a set of events {E\...En}, 
where at each event Ej a measurement of a single qubit Pauli 
operator <7j G {er , er 3 } can be made. For a particular 
choice of Pauli operators {er^ }" =1 , we take ({cr^ }" = i) to be 
the expectation value of the product of the result of these mea- 
surements. Then we can define a pseudo-density matrix 
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If the measurement events E\...E n are spacelike separated, 
then R reduces to the standard n qubit density matrix. How- 
ever, as there is no notion of separate systems inherent in the 
definition of R, it allows us to describe correlations between 
measurement events which are not spacelike separated, for ex- 
ample encapsulating the possibility of multiple measurements 
made at different points in time on a single system. This is 
a generalization of the notion of a quantum state extended 
across spacetime, rather than the usual restriction to some 
fixed time. We note that Isham iTPJl also considered intro- 
ducing a direct product structure into temporal correlations, 
but within a completely different context. 

This pseudo-density matrix inherits some properties of a 
standard density matrix. Firstly, it is Hermitian, since it is 
defined as a sum over Pauli operators, and secondly it has unit 
trace, since the expectation value ({/,...,/}) = 1. Lastly, 
the expectation value for the product of any choice of local 
measurements {cr i;j }™ =1 is given by 



<K.}? =1 )=Tr 




(5) 



As with a standard density matrix, we can also trace over sub- 
systems to produce a reduced pseudo-density matrix, defined 
only over the remaining events. 

All density matrices are positive semi-definite matrices 
with unit trace, and any matrix satisfying these requirements 
can be interpreted as a density matrix. The main difference be- 
tween R and a standard density matrix, then, is that R is not 
necessarily positive semi -definite. To see this, we consider the 
case of a single physical qubit with two separate measurement 
events. We take the qubit to be initially in the state |0) and as- 
sume that evolution between measurement events corresponds 
to the identity operator. In this case the expectation values are 
all zero, except for ({1,1}), ({X,X}), ({Y,Y}), ({Z,Z}), 
({Z, I}), and ({I, Z}), which are all equal to one. From these 
expectation values, we obtain a pseudo density matrix 
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which has eigenvalues { — ^,0, 5, 1}. The existence of nega- 
tive eigenvalues implies that R is not positive-semi definite. 

Any R which is positive semi-definite can be interpreted 
as a regular density matrix, for which it is possible to dupli- 
cate the correlations present with spacelike separated quantum 
systems. However, when R has negative eigenvalues, it can- 
not be interpreted as a regular density matrix, implying that 
the measurements cannot be spacelike separated, and hence 
there must exist a causal relation between events. 

The causal relationship embodied in certain pseudo-density 
matrices has may similarities to another form of uniquely 
quantum correlation, entanglement, and we can define an 
analogous measure of causal correlations. In order for a func- 
tion f(R) to be considered a causality monotone we require 
the following criteria to hold: 

1- f(R) > 0, with f(R) = if R is completely positive, 
and /(-R2) = 1 for any R2 obtained from two consecu- 
tive measurements on a single qubit closed system, 

2. f(R) is invariant under unitary operations, 

3. f(R) is non-increasing under local operations, and 

4. EiPif(Ri)>f(£iPiRi)- 

These criteria correspond almost exactly to the criteria for an 
entanglement monotone lfT4l[T5l . except that criterion three is 
somewhat weakened. An entanglement monotone is required 
not to increase on average under local operations and classi- 
cal communication, however any processing based on classi- 
cal communication would constitute a causal relationship and 
hence is excluded. 

As we have shown, any pseudo-density matrix which em- 
bodies some form of causal relationship must have at least one 
negative eigenvalue. Since such matrices are Hermitian (and 
hence have real eigenvalues) and have unit trace, it follows 
that the trace norm is strictly greater than one. On the other 
hand, if all eigenvalues are positive, the trace norm is exactly 
one. 

This leads us to define a measure based on the trace norm, 
ftr(R) = ||-R||tr — 1- As we have seen, ||i?||t r > 1 for all 
valid pseudo-density matrices, and hence ftr(R) > 0. Fur- 
ther, ftr(R) = trivially for all positive semi-definite R, and 
from the previous example it is clear that ftriRi) — 1 for 
at least one choice of i? 2 . Since the trace norm is unitarily 
invariant, the first and second criteria for f tr to be a causal- 
ity monotone are satisfied. Similarly, by applying Stinespring 
dilation to represent local quantum operations as unitary op- 
erations on a larger Hilbert space, the third criterion follows 
directly since the trace norm is non-increasing under partial 
trace. The final criterion follows from the triangle inequal- 
ity since || Y Jl Pi R i\\tr - 1 < Y,iPi(\\ R t\\tr ~ 1) and hence 
/*r(X)j Pi-Ri) < Y.iPiftr{Ri)- Thus f tr is a causality mono- 
tone. 

Naively it would appear that non-destructive single qubit 
measurements are necessary in order to perform the multi- 
event measurements required for tomography of a pseudo- 
density matrix. This would rule out the possibility of re- 
constructing a pseudo-density matrix in either NMR or quan- 
tum optics, two of the most established testbeds for quantum 
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FIG. 1: Eigenvalues of R and value of /„ as a function of T wa i,. In (a) the system starts in a pseudo-pure state and undergoes dephasing noise, 
while in (b) the system starts in a mixed state and undergoes depolarising noise. The circles indicate data points obtained from experiment, 
while the solid lines indicate the best fit for the relevant theoretical models. These models each take 3 parameters to describe the initial state of 
the system and either 1 and 3 parameters, respectively, to parametrize the noise. The red region indicates the time period in which all resulting 
pseudo-density matrices are acausal. 



physics. Fortunately, however, it is possible to circumvent 
the limitations imposed by ensemble measurements by mak- 
ing use of an ancilla qubit to record the parity of the local 
Pauli measurement results. Thus it is possible to recover their 
product by measuring a single spin, similar to the approach 
advocated in |fl6l . 

The simplest system for which R can have negative eigen- 
values contains two measurement events. These can be made 
either on the same qubit or seperate qubits. However in order 
to observe both causal and acausal correlations in the current 
generation of experiments we focus on measurements sepa- 
rated by a variable time on a single qubit, as in the circuit 
below which accomplishes the measurement {<7i, 02}. 
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Here U ai{2) is the unitary operation mapping the ±1 eigen- 
state of the Pauli operator cr 1 ( 2 ) onto the ±1 eigenstate of 
Z. Between measurements we allow a period of free evolu- 
tion, during which the primary qubit undergoes decoherence, 
and we calculate a pseudo-density matrix Rt wA for a range of 
waiting times. 

NMR experiments were performed on a Varian Unity Inova 
spectrometer with a nominal 1 H frequency of 600 MHz using 
a HF{CP} probe with pulsed field gradients. The NMR sam- 
ple comprised 13 C-labelled sodium formate dissolved in D2O 
at 20°C, providing a heteronuclear two-spin system. The : H 
spin was used as the primary qubit and the 13 C spin as the an- 
cilla. Both spins were placed on resonance, so that the Hamil- 
tonian took the form of a spin-spin ZZ coupling of 194.7 Hz, 
and the Bi field strengths were adjusted to give nutation rates 
of 12.5kHz. The measured relaxation times were Tj=7.8s 
and T 2 =3.2 s for X H and Ti=16.3 s and T 2 =6.7 s for 13 C. An 



inter-scan delay of 60s ensured that the spin system began 
each experiment close to its thermal state. 

Quantum logic gates were implemented using standard ap- 
proaches ifTTl rrSl . Single qubit rotations in the XF-plane 
were implemented using BB1 composite rotations lfT9l l20l . 
while Z-rotations were implemented as frame rotations ETI 
which were propagated through the pulse sequence E2l to 
points where they could be dropped. Pseudo-pure two-qubit 
states were prepared using the method of Kawamura et al. 
||23l ; for pseudo-pure single qubit states the thermal state was 
used directly. NMR spectra were processed using home writ- 
ten software and the intensity of the 13 C doublet determined 
by combining separate integrals for the two components; all 
integrals were normalised using a reference spectrum. 

Instead of using natural decoherence during T wa ; t , control- 
lable dephasing of the primary qubit was implemented using 
the diffusive suppression of pulse field gradient spin echoes 
E4l as described by Cory et al. ||25l . This can be converted 
to controlled depolarization by using single qubit rotations to 
apply the dephasing around the X, Y and Z-axes in turn. 
This process also dephases the ancilla qubit, but leaves its Z- 
component unaffected as the ancilla does not experience the 
single qubit gates. 

Figure [T] shows the results of our NMR experiments, plot- 
ting the eigenvalues of the pseudo-density matrices as a func- 
tion of time, along with the corresponding f tr , in each of two 
settings. Figure [T}\ shows the results of purely dephasing 
noise acting on an initial state pseudo-pure state |0). Here 
the pseudo-density matrix starts with a single negative eigen- 
value, which tends towards zero from below as the waiting 
time is increased. The pseudo-density matrix never becomes 
positive semi-definite (and hence acausal) because the deco- 
herence brings it towards a matrix which is rank deficient, 
and so the minimum eigenvalue approaches, but never quite 
reaches, zero. 
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In order to observe a sharp transition between causal and 
acausal pseudo-density matrices it is necessary both to start 
with a mixed initial state, and to allow depolarizing deco- 
herence, which is the case considered in Figure [Tj3. Now 
we observe a transition between causal and acausal pseudo- 
density matrices as the minimum eigenvalue crosses the zero 
threshold, a phenomenon reminiscent of entanglement sudden 



death |26l . 
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